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We report the first measurements of the dynamics of liquid germanium (l-Ge) by quasi-elastic 
neutron scattering on time-of-fiight and triple-axis spectrometers. These results are compared with 
simulation data of the structure and dynamics of l-Ge which have been obtained with ab initio density 
functional theory methods. The simulations accurately reproduce previous results from elastic 
and inelastic scattering experiments, as well as the q-dependence of the width of the quasi-elastic 
signal of the new experimental data. In order to understand some special features of the structure 
of the liquid we have also simulated amorphous Go. Overall we find that the atomistic model 
represents accurately the average structure of real l-Ge as well as the time dependent structural 
fluctuations. The new quasi-elastic neutron scattering data allows us to investigate to what extent 
simple theoretical models can be used to describe diffusion in l-Ge. 
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I. INTRODUCTION 

The investigation of the structure of the condensed 
phases of germanium is of importance for fundamental 
science since this material exhibits polymorphism in the 
solid state and pronounced changes in density and bond- 
ing upon melting. Research in this field is also relevant 
to several applied physics problems and semiconductor 
technology. 

In its solid phase germanium shows a large degree of 
polymorphism. In its stable crystalline phase at ambient 
pressure, the crystal has a diamond structure in which 
each Ge atom is surrounded by four covalcntly bonded 
first neighbors in a tetrahcdral formation. This phase of 
Ge has a very low density with respect to a close-packed 
structure and is semiconducting. Upon application of 
pressure, the tetrahedrally bonded network is disrupted 
and the number of neighbors and the density increase. 
Around 100 kbar, a phase transition to the metallic (3- 
tin (or white-tin) structure occurs and at even higher 
pressures an hexagonal phase and a close-packed phase 
are founds. 

At standard pressure the melting temperature of crys- 
talline Ge is T m = 1210.4 K. At the melting transition 
the tetrahcdral network is disrupted and the average co- 



ordination number increases from 4 to about 7, which is 
still low compared to other liquid metals which typically 
have coordination number between 9 and 12. Liquid ger- 
manium is metallic^. Germanium can also be produced 
in an amorphous phase in which distorted tetrahedra are 
connected in a continuous random network with defects. 
Amorphous germanium is semiconducting. 

Experimental investigations of the density and bond- 
ing changes in different phases of germanium are numer- 
ous, using both X-ray^i£ and neutron diffraction&i&SJJii 
techniques. The static structure factor shows features 
that differ from those found in simple liquids: The first 
diffraction maximum is unusually low and has a shoulder 
on its high momentum transfer side. Thus a hard-sphere 
model does not seem to be appropriate for describing 
the arrangement of atoms in l-Ge and several models 
have been proposed based on the short-range order found 
in white-tin (also called (3— tin)^, on the coexistence of 
randomly distributed, molten metal-like atomic arrange- 
ment and covalent crystal-like atomic arrangement^, on 
a mixture of fourfold and highly coordinated metallic 
arrangements^, or on a quasicrystallinc model of i-Ge£. 
However none of these accounts for all of the structural 
features of l-Gc. 

In order to investigate the dynamics of l-Ge, Hosokawa 
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and co- workers 11 used inelastic X-ray scattering to mea- 
sure the spectra for momentum transfers q from 0.2 
to 2.8 A -1 . The phonon dispersion curve extracted 
from the low q spectra matches the hydrodynamic sound 
velocity^. The study of the relationship between S(q) 
and T q , the half width of S(q, lo) at half maximum, shows 
that at higher q the measured spectral width cannot be 
described by the relationship given by de Genncs for sim- 
ple dense fluids 1 ^ nor the expression for spectra of dense 
hard-sphere fluids given by Cohen and co-workers^ 4 -. Ac- 
cording to Hosokawa et al~^, this discrepancy supports 
the view that neither the structure nor the dynamics of l- 
Ge can be described using a single pair interaction. How- 
ever, Ashcroft^ proposed a cluster model with transient 
covalent structures in the liquid which seems to be con- 
sistent with the experimental results. In this model, the 
first maximum in S(q) reflects the inter-cluster correla- 
tion while the shoulder on the high q side and the re- 
maining oscillations are associated with the atom-atom 
contributions. Regarding the experimental results on the 
dynamics, Hosokawa et al. put forward the idea that the 
narrowing of the quasi-elastic line at the position of the 
maximum in S(q), followed by an increase of T q at q val- 
ues that correspond to the location of the shoulder in 
S(q), indicates that the covalent structures are diffusing 
slowly while individual atoms are subject to rapid trans- 
lational motion. 

Also on the theoretical side, much effort has been de- 
voted to modeling and reproducing both structural and 
dynamic features of l-Ge. In 1985, Stillingcr and We- 
ber proposed a model potential function comprising both 
two- and three-atom contributions to describe the inter- 
actions in solid and liquid states of Sii£. Molecular dy- 
namics simulations of l-Si using this potential give struc- 
tural properties in good agreement with the measure- 
ments. Subsequently this potential was also used for sim- 
ulating amorphous germanium^. Later on the potential 
devised for l-Si was used with a different set of param- 
eters which were obtained by fitting crystalline as well 
as liquid Ge phases 1 ^, the corresponding static structure 
factors and pair correlation functions therefore being in 
good agreement with the experimental results. However, 
it was found that the computed diffusion coefficient is 
larger than the experimental values by a factor of 2. 

There have also been a number of ab initio molecu- 
lar dynamics (AIMD) simulations of Z-Ge, starting in 
1993 with the work of Kresse and Hafner who found 
a quite good agreement between simulated and experi- 
mental pair correlation functions and electronic densities 
of states 1 ^. Subsequently they extended their study to 
the liquid metal-amorphous semiconductor transition^. 
Their results for the pair correlation function and the 
static structure factor are in good agreement with ex- 
perimental measurements for the amorphous phase. Ac- 
cording to these authors, the results for the liquid in- 
dicate a broad homogeneous distribution of local bond- 
ing configurations and show that both classes of mod- 
els proposed earlier (mixture of fourfold coordinated and 



highly coordinated metallic arrangements^ and both six- 
fold coordinated and metallic arrangements^) are unre- 
alistic. The simulation of supercooled and amorphous 
Ge shows the increase of the local tetrahedral order with 
decreasing temperature. The self-diffusion coefficient ex- 
tracted from the simulated trajectory close to the melt- 
ing point, D = 1.0- 10 -4 cm 2 /s, is in good agreement 
with the measured value. Similar results about the struc- 
ture and diffusion of l-Ge were found by Takeuchi and 
Garzor*2i using the ab initio Car-Parrincllo molecular- 
dynamics scheme^. A further study was performed by 
Godlevsky et al*^ which yielded similar results for both 
the structure factor and the pair correlation function. 
Kulkarni and collaborators performed similar simulations 
and studied the influence of the temperature between 
1250 K and 2000 K on the structural properties of the 
liquid^ 4 -. These simulations show that liquid Ge is a good 
metal but has some special short-range order arising from 
residual covalent bonding. Only at the highest temper- 
atures, l-Ge seems to evolve toward a more conventional 
close-packed liquid metal. Calculations of the dynamic 
structure factor from an ab initio MD simulation were 
performed by Chai and co-workers2£. They show good 
agreement with the experimental data, although the sta- 
tistical noise due to the small number of atoms (64) is 
rather large and has to be convoluted with a Gaussian of 
width 2.5 meV representing the experimental resolution 
function. 

As can be seen from these experimental and theoretical 
studies, experimental investigations of l-Ge are presently 
not able to determine the microscopic, atomic arrange- 
ments in the liquid, but the experiments provide the data 
which is needed to validate theoretical models and sim- 
ulations. In this paper, we present a new study of liquid 
germanium using both quasi-elastic neutron scattering 
(QENS) experiments and AIMD simulations. The QENS 
technique provides much clearer insight into diffusive dy- 
namics than inelastic X-ray scattering in an intermediate 
q range due to the higher resolution and approximately 
Gaussian lincshapc of the spectrometer resolution func- 
tion. The QENS signal is typically most prominent at the 
g-rangc around the first sharp diffraction peak. QENS 
data therefore relates directly to local structural fluctua- 
tions on a length scale which is most relevant to the on- 
going discussion about the structure of l-Gc. Section HTl 
presents details and results of the QENS measurements. 
Section IIIII is dedicated to the ab initio simulations of 
the structure and dynamics of l-Ge and to the compar- 
ison with new and existing experimental results and in 
Section lTvl we summarize and conclude the present work. 



II. NEUTRON SCATTERING STUDY OF l-Ge 

QENS measurements on l-Ge were performed on cold 
neutron, time-of-flight (Focus at the Swiss spallation neu- 
tron source SINQ, PSI, Switzerland2&2!i2£) and triple- 
axis (4F1 at Laboratoire Leon Brillouin, France^) spec- 



3 



tromctcrs. Note that due to the high sound velocity, 
c s = 2682 m/si£, the acoustic phonon of l-Ge could not 
be measured. 



A. Sample preparation and environment 

Solid Ge samples were obtained either as a powder (ob- 
tained commercially with a purity of 99.99% in weight) or 
as pure solid pieces (less sensitive to oxidation and higher 
packing factor). Solid Ge was inserted into a cylindrical 
cell (6 mm in diameter) made of vitreous silica which was 
then sealed under vacuum. Sample preparation was per- 
formed either under nitrogen atmosphere or vacuum, in 
order to avoid contamination of dioxygen (which would 
lead to formation of GeO or GeC>2 at high temperatures). 
The silica cell was suspended from the end of the sample 
stick in a niobium envelope. The melting temperature 
of l-Ge is T m = 1210.4 K. During the neutron scattering 
experiments, Ge was heated in standard furnaces. Af- 
ter melting, the height of the sample was between 2 and 
5 cm, depending on the initial solid form of Ge. 



B. Experimental setup and data analysis 

In both experiments, the sample temperature was com- 
prised between 1260 K and 1400 K, with a sample of 
6 mm in diameter. 

On Focus, the measurements were performed at inci- 
dent wavelengths of 3 A and 4 A. Momentum transfer 
ranges from 0.35 to 3.8 A -1 at A; = 3 A. Background 
measurements with comparable statistics were performed 
at both wavelengths for the furnace, empty sample cell 
and niobium envelope. The instrument resolution func- 
tion was determined using a cylindrical vanadium sam- 
ple of diameter 10 mm. The clastic energy resolution 
at incident wavelength of 3 A and 4 A was of the or- 
der of 500 ^eV and 200 fxeV, respectively. The data 
were processed with standard software at the instrument 
(NATHAN). 

On 4F1 constant q scans were performed between 
q = 0.3 A" 1 and q = 3.9 A -1 , although the QENS signal 
could only be measured unambiguously for q > 1.5 A -1 . 
The data were processed with the standard 4F1 program. 



C. Experimental results 

The background signal was almost as strong as the sig- 
nal from l-Gc. A slight misalignment between the sam- 
ple and empty cell measurements meant that the data 
around the elastic peak had to be discarded. The dif- 
ference spectra were fitted using a single Lorentzian for 
each q (see Fig. [T]for Focus results). The QENS inten- 
sity increases and the width decreases as q approaches 
2.5 A -1 . The triple- axis measurements give similar re- 
sults (see Fig. [J). The evolution of the QENS signal 
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FIG. 1: Time-of-flight measurements of S(q,ui) for l-Ge at 
different q values (\i = 3 A; T = 1260 K); different spectra 
are shifted vertically by 20 counts. The solid lines show the 
fitted Lorentzians. 
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FIG. 2: Triple-axis measurements of S(q, u>) for l-Ge at dif- 
ferent q values around the position of the first maximum in 
S{q) (k f = 2.662 A" 1 , T = 1300 K). 



with q will be discussed below with the results from the 
simulation. 



III. AB INITIO SIMULATION STUDY OF l-Ge 
A. Simulation procedure 

All calculations were performed with the Vienna Ab 
initio Simulation Package (VASP) 30 i 31 , which is based 
on density functional theory. We carried out most sim- 
ulations of l-Ge for an ensemble of 64 atoms in a cubic 
supercell with periodic boundary conditions. The simu- 
lated density was p = 5.53 g/cm 3 , which corresponds to 
a box of length 11.175 A. 

The atomic dynamics was described in the micro- 
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canonical (N,V,E) ensemble. Ultrasoft Vandcrbilt 
pseudo-potentials^ and the generalized gradient ap- 
proximation (GGA) with the functional PW91 for the 
exchange-correlation energy were used. The Brillouin 
zone was sampled by a 3 x 3 x 3 Monkhorst-Pack mesh of 
fc-points^ 3 -. The wave functions were expanded in a ba- 
sis of more than 32000 plane waves, corresponding to an 
energy cut-off of 139.2 eV. One MD run was performed 
in a cell containing 200 atoms for which the fc-point grid 
contained only the gamma point. 

The starting configurations for the ab initio simula- 
tions were generated from classical MD simulations using 
the Stillingcr and Weber potentiate. For this purpose, 
systems of 64 Ge atoms were equilibrated at the desired 
temperatures (three independent samples were extracted 
from a 50 ps run at high temperature (1700 K); these 
samples were then equilibrated at the desired tempera- 
ture by velocity scaling during 10 ps). 

After switching to the ab initio (N,V,E) MD, and an 
equilibration time of 5 ps, three runs of 30 ps using a 
time step of 3 fs were performed at three different tem- 
peratures corresponding to the liquid (1390 K, 1080 K) 
and supercooled (680 K) states, all at the same density. 
A similar simulation scheme was used to perform a liq- 
uid MD run at 1390 K for the cell containing 200 atoms, 
the NVE trajectory spanning 15 ps. Three independent 
configurations of the trajectory at 680 K were quenched 
at a cooling rate of 68 K/ps in order to obtain three 
amorphous configurations of the system. 



B. Numerical results 

The mean-squared displacement (MSD) of atoms is one 
of the dynamic quantities that allows to check whether 
the system is a liquid or not. The MSD is computed from 
the relation 



dr 2 (t) = ([r l (t)-r l (0)} 2 ) 



(1) 



and allows the diffusion constant D to be determined via 
the Einstein relation : 



D = lim 



! (t) 



6f 



(2) 



In Fig. [3] we compare the values for D as obtained from 
the simulations with experimental data from Pavlov and 
Dobrokhotov 34 , For T around 1400 K, very similar val- 
ues of D are obtained from experiment and ab initio MD, 
while according to Ref. [H classical MD overestimates D 
by a factor of 2, which indicates that from this point 
of view the ab initio calculation is indeed more accu- 
rate than simulations with classical force fields. Having 
performed liquid simulations at 3 temperatures, we are 
able to observe an Arrhcnius behavior of the diffusion 
constant with an activation energy around 0.72 eV. This 
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FIG. 3: Arrhenius plot of calculated and experimental diffu- 
sion constant D (experimental data from Ref. [Hj). Note that 
the vertical scale is logarithmic. 



quantity corresponds to the average energy needed by a 
Ge atom to escape from its local environment and un- 
dergo diffusion and it is characteristic of a strong liquid. 
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FIG. 4: Static structure factor of l-Ge: Simulated data for 
T = 1390 K (dashed line), T = 1080 K (dotted line), and 
T — 680 K (full line) and experimental data by Salmon^ at 
T = 1270 K (triangles). 

The structure of the liquid can be characterized by the 
static structure factor which is defined as 3 ^: 



N N 



(3) 



3 = 1 1=1 



It was computed using, at each temperature, 10000 atom- 
istic configurations. The average is also performed over 
all q vectors of the same magnitude. The results are 
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FIG. 5: Pair correlation function of l-Ge: Simulated data: 
T = 1390 K (dashed line); T = 1080 K (dotted line) and 
T = 680 K (full line) and neutron scattering data by Salmon 10 
at T = 1270 K (triangles). 

shown in Fig. [4] together with the experimental result by 
SalmoniS. 

At T = 1390 K the main sharp peak occurs at 
q = 2.53 A -1 which is very close to the experimental 
values (q m ax ~ 2.50 A -1 )££*!£. At lower temperatures, 
the peak is shifted to slightly lower q values, which is 
consistent with measurements in the temperature range 
from 1270 K to 1820 K2. The characteristic feature of l- 
Ge shown by neutron diffraction studies^, the distinct 
shoulder at q = 3.27 A" 1 , i.e. on the high-g side of 
the first sharp peak, is reproduced well by the simula- 
tion. As T decreases, the height of the main peak in- 
creases and the shoulder becomes more distinct. These 
two features were described at higher temperatures in 
the simulation work by Kulkarni et al. who found that 
at T = 2000 K the shoulder disappears^. The second 
peak of the static structure factor occurs at 5 A -1 which 
is slightly shifted from the experimental value reported 
by Salmon^ who found q = 5.1 A -1 (see also Waseda£), 
although the experimental value seems to have an uncer- 
tainty of « 0.1 A" 1 !^. 

The pair correlation function g(r), defined &s~ 

N N 
i=l j^i 

allows the structure of the system in real space to be 
characterised. A comparison of simulated and experi- 
mental g(r) is given in Fig. [5] In particular we see that 
g(r) hardly shows a second peak at distances beyond the 
first peak, in marked contrast to the situation found in 
simple liquids such as hard spheres. The structure of l- 
Ge is significantly different from that of simple liquids. 
In addition the height of the main peak increases with 
decreasing temperature but the structure at larger dis- 
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FIG. 6: S(q) for supercooled and amorphous Ge: Comparison 
between simulation (supercooled state at T = 680 K (open cir- 
cles) and amorphous state (open diamonds)) and experimen- 
tal data, from Ref. |37| . as obtained from neutron scattering 
with an amorphous sample (filled diamonds). 



tances remains basically unchanged. Indeed both S(q) 
and g(r) show that structural changes over a wide range 
of temperatures from 680 to 1390 K are generally rather 
small. 

Structural changes become more significant if the liq- 
uid phase is quenched to produce an amorphous phase. 
This can be seen in Fig. [5] in which we show S(q) for 
the supercooled liquid at 680 K and the amorphous 
state. The amorphous phase was obtained by cooling 
down, over 10 ps and at constant volume, three indepen- 
dent configurations of the supercooled state at 680 K to 
T = K. Thus for this state S(q) was calculated by aver- 
aging over only three configurations, which is the reason 
for the substantial noise in the simulation data. In S(q), 
the first peak is shifted to higher q reflecting the high, 
liquid density that we have imposed on the amorphous 
structure. At higher g, the agreement between experi- 
ment and simulation is much better. Looking at snap- 
shots of the system in the amorphous state, one clearly 
sees a tetrahedral structure in which a Ge atom occupies 
the center of a tetrahedron. The calculated g(r) for the 
amorphous state (not shown here) shows a main peak 
shifted to lower interatomic distances (as compared to 
the liquid state) while the intermediate peak near 4 A be- 
comes much more pronounced, in agreement with experi- 
mental findings 3 -^. The main peak in g(r) corresponds to 
the interatomic distance between the center and a vertex 
of the tetrahedron while the higher interatomic distance 
at 4 A is associated with the distance between two ver- 
tices of the tetrahedron (see also Fig. [7]). Despite im- 
posing the higher liquid density on our amorphous struc- 
tures, the shortest Ge-Ge distances are overestimated by 
several 0.1 A. 

Given the nature of the many-body interactions in ger- 
manium and the resulting open structure, it is interesting 



FIG. 7: (Color online) Calculated structure of amorphous Ge at T = K. Left: Snapshot of the whole simulation box (with 
small spheres representing the corners of the box). Middle: Interatomic distances. Right: Bond angles. 
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FIG. 8: Angular distribution (left) and coordination number 
(right) for liquid (top), supercooled (middle) and amorphous 
(bottom) phases of Ge for r m in the vicinity of the first min- 
imum of g(r). Circles: r m = 2.8 A; Squares: r m = 3.0 A; 
Diamonds: r m = 3.2 A; Triangles r m — 3.4 A. The quanti- 
ties for the liquid, supercooled, and amorphous phases are, 
respectively, averaged over 9000, 6300 and 3 configurations. 



to compute the angular distribution of atoms. For this 
we define P(9,r m ) as the normalized distribution of an- 
gles defined by vectors going from a reference atom i to 
two neighbors of i which are within a radius r m of the 
reference atom. In Fig. [8] we show P(9,r m ) for several 
values of r rn around the minimum in g(r). In the liquid 
this minimum is not very pronounced, see Fig. [51 leading 
to an uncertainty regarding the definition of the near- 
est neighbor shell. In contrast, in the amorphous state, 
the first minimum in g(r) (calculated, not shown here) 
is clearly defined allowing an unambiguous definition of 
nearest neighbors. We note that any bond angle distri- 
bution plot will tend to have a maximum close to 90 ° and 
that changes in the distribution as a function of radius 
and temperature are therefore particularly meaningful. 

For liquid and supercooled systems, we observe a broad 
peak around 6 = 98 ° if r m = 2.8 A. With increasing 
r m , this peak tends to shift to smaller angles, around 
= 90 ° . For the amorphous phase, the corresponding 
peak is observed at slightly larger angles, around 105 ° for 
r m = 2.8 A, indicative of a slightly more regular tetrahe- 
dral environment in the amorphous phase, even at higher 
pressure as imposed by the liquid density. With increas- 
ing r m , a new peak appears for angles between 55 ° and 
60 " in the liquid and supercooled phases. At T — 680 K, 
this feature is more pronounced than at higher temper- 
atures. In the amorphous phase, this peak appears at 
smaller angles, about 50 ° . A bond angle of 60 ° is char- 
acteristic of metallic bonding, since this value represents 
the locally most closely packed structure of the system. 
The broad peak around 100 ° is typical of flattened tetra- 
hedra, the symmetric tetrahedral angle being 109.5 ° , see 
also Fig. [3 

Regarding the coordination number, most atoms have 
either 4 or 5 nearest neighbors in the amorphous phase 
(see right panels of Fig. [5]). On the other hand, in the 
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liquid and supercooled phases, for typical values of the 
first minimum in g(r), 3.4 A, the coordination number is 
between 6 and 7. 

We note that obviously the coordination number in- 
creases with increasing radius of the coordination sphere 
for both liquid and supercooled systems, which is consis- 
tent with the so-called (3-t'm hypothesis^. This hypothe- 
sis states that in /-Ge, each atom has four nearest neigh- 
bors on a flattened tetrahedron and two other neighbors 
slightly further away on the normal to the plane of the 
flattened tetrahedron (see Fig. [9] for a schematic repre- 
sentation of the local environment of j3— tin). Denoting 
by ro the main interatomic distance reported from g(r), 
i.e. ro ~ 2.5 A, the /3-tin configuration would lead to in- 
teratomic distances of the order of ro (between the center 
and a vertex of the tetrahedron), v^ro ~ 3.53 A (be- 
tween two adjacent vertices of the tetrahedron), and 
2ro ~ 5 A (between two opposite vertices of the tetra- 
hedron) together with distances slightly higher than ro 
associated with the two further neighbors. 

Snapshots of the simulated structure of amorphous Ge 
are shown in Fig. [7] On the left, the 64 atoms are rep- 
resented together with the corners of the simulation box. 
The bonds between atoms correspond to two covalent 
radii (1.22 A for Ge) plus a bond tolerance of 0.4 A so 
a bond is drawn for interatomic distances smaller than 
2.84 A. The two other plots show a zoom into this struc- 
ture. Three "tetrahedra" are represented together with 
some interatomic distances and bond angles. The cova- 
lent bond length of Ge corresponds to the main peak in 
g{r) around 2.5 A. On the other hand, the distance be- 
tween two neighboring vertices of a tetrahedron ranges 
from 3.5 A to 4.5 A, which accounts for the broad peak 
around 4 A in g(r) (calculated, not shown here). Looking 
at the bond angles, their values are mainly distributed 
between 85 ° and 130 ° with an average value of 105 ° , 
close to the tetrahedral value (see Fig. [SJ. At higher 
interatomic distances (see distribution for r m = 3.4 A in 
Fig. , another peak appears around 50 ° . This one 
is due to the angles formed by the two following vectors: 
The first one goes from a vertex to the center of the tetra- 
hedron while the second one goes from the same initial 
vertex to another vertex of the same tetrahedron. This 
is also shown in Fig. [7] 

Having established that the liquid simulations do in- 
deed describe well the corresponding experimental sys- 
tem and having studied the structural aspects, we now 
investigate the dynamic features of the system. For this 
we start with the intermediate scattering function which 
is defined as^ 

JV N 

%^-^^(cxp[, 9 .(r J (()-r t (0))]) . (5) 

3=1 k=l 

Since the system is isotropic, we have performed an av- 
erage over the different orientations of q. 

Figure [TO] shows the (normalized) F(q, t) of /-Ge at 
1390 K and for low q values from the large simulation 



cell. The oscillations are associated with the presence 
of the acoustic mode at these wave-vectors. The pe- 
riod of these oscillations decreases with increasing q and 
also they become more damped, which is directly associ- 
ated with the shape of the dispersion and damping of the 
acoustic phonon as a function of q, as will be discussed 
in more detail below. 

The dynamic structure factor S(q,u>) was computed 
by a Fourier transform in time of F(q, t)^. In a periodic 
system of cell side L only the reciprocal lattice points 
(n x ,n y ,n z ).2iT/L are accessible. The first q-point corre- 
sponds to a spatial correlation over a distance L, which 
is not meaningful in a cell of side L. g-points bigger than 
An / L correspond to correlations over distances shorter 
than L/2, which are uniquely defined. In the interme- 
diate range there is a progressive increase in meaning- 
ful correlations which correspond to distances between 
atoms close to the centre and those close to the apexes 
of the cell. The lowest unambiguous q-values, calcu- 
lated in 0.1 A -1 ^-strips, for the 200 atom cell start at 
0.7 A -1 and at 1.1 A" 1 for the 64 atom cell. At smaller 
g, <7 < 1 A -1 , S(q, u>) shows a well defined acoustic mode, 
the position of which can be compared directly with the 
corresponding experimental results from the inelastic X- 
ray scattering of Ref. ITU Such a comparison is made 
in Fig. [11] and it shows a very good agreement between 
experiment and simulations. The width of the acous- 
tic phonons increases with frequency, both in the exper- 
iment and in our simulations. At 1.2 A" 1 the phonons 
are hardly visible in the experimental spectrum and in 
the simulations, the width is of the same order as the 
frequency. Accordingly we estimate the error in the fre- 
quency to be a few meV. We also mention that the disper- 
sion of the acoustic mode does not show any significant 
variation with T in the temperature range investigated. 

At higher q values, q > 1 A -1 , the QENS signal domi- 
nates S(q,u)) (see Fig. fl2]> . Fitting the quasi-elastic line 
of S(q, uj) with a Lorentzian allows to determine its half- 
width at half maximum (HWHM), the (/-dependence of 
which is shown in Fig. [13] together with the correspond- 
ing structure factor. Unconstrained fits of the Focus data 
lead to an underestimation of the spectral width at low q 
(< 1.8 A -1 ), where the QENS intensity is weak, and also 
at higher q (> 3.2 A -1 ), where only the negative energy 
transfer part of the spectrum is present, and therefore 
we show the Focus data only for 1.8 A -1 < q < 3.2 A -1 . 
The simulated widths are obtained from the simulation 
at 1390 K by fitting S(q, to) with a Lorentzian. From 
Fig. [13] we recognize that the data from the simulation 
agrees very well with the one from the experiment. In 
particular we see that in the simulation as well as in the 
experiments a narrowing of the quasi-elastic peak is ob- 
served for q = 2.4 A -1 and between q = 2.8 A -1 and 
3.2 A -1 , which corresponds to the main peak and the 
shoulder in the static structure factor, respectively (see 
line in Fig. [T3"|) . This behavior can be interpreted as de 
Gennes narrowing^, i.e. the fact that if one considers 
the dynamics on length-scales that are close to the one 




FIG. 9: Schematic representation of the /3-tin configuration. Atom 1 is at the center of the flattened tetrahedron formed 
by atoms 2, 3, 4, and 5. Atoms 6 and 7 are the second nearest neighbors. Distances are: di2 = di3 = di4 = dis « ro ; 
d 2 3 = d 3 A = d 45 = d 52 ~ V2r ; d 2 4 = ^35 ~ 2r . 
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FIG. 10: Collective dynamic properties of l-Ge at low q val- 
ues: Simulated intermediate scattering function F(q, t) at 
T = 1390 K normalized by its value at t = 0. 

of the local structure of the liquid, the structure is basi- 
cally preserved and thus there is only a relatively small 
loss of memory and consequently a small broadening of 
the dynamic structure factor. 

The g-dependence of the width of the quasi-elastic 
peak can be compared with theoretical predictions based 
on different models. The simplest expression for the 
HWHM for a dense monoatomic liquid is that the half- 
width ojq is the square root of the normalized second 
frequency moment^: 

H 2 ) = ^T? 2 ■ ( 6 ) 

In a somewhat refined approach. Cohen et alr^ pre- 
dicted for a dense hard-sphere fluid the half-width u>h to 
be given by 



q I i I i I i I i I i I 

0.25 0.5 0.75 1 1.25 
q (A"') 

FIG. 11: Dispersion of the frequency of the acoustic phonon as 
a function of q: Comparison between inelastic X-ray scatter- 
ing data at T = 1250 K from Ref.[ll] (squares) and simulation 
at T = 1390 K (circles for the 200 atom cell, triangles for the 
64 atom cell). 



wh(q) = "57^"%)' ( 7 ) 

where 

d(q) = [l-j (qa) + 2j 2 (qa)}- 1 . (8) 

Here De is the self-diffusion coefficient of the hard-sphere 
fluid within the Enskog theory, a is the diameter of 
the hard spheres, and jo(x) and ]2{x) are the zeroth- 
and second-order spherical Bessel functions, respectively. 
Note that equation {7} is only valid for intermediate val- 
ues of q, 1 < qa < a /I, where a /I S> 1 with I the mean 
free path of a particle which, since we arc considering a 
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FIG. 12: Simulated dynamic structure factor for /-Ge at 
1390 K for different q values (Symbols: Simulation; Line: Fit 
with a Lorentzian). 




q(A ) 



FIG. 13: Half width at half maximum of S(q, u>) as a function 
of q together with the static structure factor. Circles: Triple- 
axis measurements (4F1); Triangles: Time-of-flight measure- 
ments (Focus); Squares: Simulation; Stars and straight line: 
Simulated S(q) at 1390 K. 



dense system, is much smaller than a. Approximating 
the Ge liquid by a system of hard spheres, we have cal- 
culated the quantity 0Jh(q) by using the static structure 
factor from the AIMD simulations at 680 K and 1390 K. 
For the size of the hard spheres, the location of the first 
peak in the pair correlation function gives, see Fig. [5l 
a = 2.8 A at T = 1390 K and a = 2.7 A for T = 680 K. 
For De, we used the value of the self-diffusion constant 
D calculated from the simulation (its dependence on the 
mass and the density can be expected to be very sim- 
ilar to that of the self-diffusion constant of the Enskog 
theory). 

These models are compared with the width of the 
quasi-elastic signal obtained from AIMD in Fig. [H] Both 
models show de Gennes narrowing at q = 2.3 A" 1 . For 
the dense monoatomic fluid model, the q-dependence of 




q(A) 



FIG. 14: Half width at half maximum of the dynamic struc- 
ture factor as a function of q for T = 1390 K. Squares: Sim- 
ulation; Diamonds: Second frequency moment ui q of a mono- 
matomic liquid, see Eq. JB); Circles: uj^ as calculated from 
the theory of dense hard-sphere fluids, see Eq. ([7]). 



uj q is reasonably good but the value of u> q is generally 
overestimated by a factor of 2. Such a systematic dis- 
crepancy has also been found for a variety of simple liq- 
uids (Ar, Kr, Lennard-Jones)^ and thus it is not surpris- 
ing that in the present system, which has a much more 
complex structure, we find a similar deviation. In con- 
trast to this the width loh predicted by the hard sphere 
model is in better agreement close to q — 2.3 A -1 but 
the increase in width away from this g-value is over- 
estimated. Also this result is in qualitative agreement 
with the ones reported by Cohen et al. for simple liq- 
uids^ which thus gives evidence that Eq. [7| is not valid 
for simple liquids and for network forming systems. 



IV. DISCUSSION AND CONCLUSION 

Liquid Ge is a system of fundamental interest for its 
structure and fluctuations. Accordingly a number of 
models have been proposed to describe the local struc- 
ture of this liquid. Due to the complexity of the system, 
models can only be validated by comparison with experi- 
mental data but it is found that static structure measure- 
ments presently available are not sufficient to distinguish 
between various plausible models. Since the local struc- 
ture of a liquid is time dependent, a more stringent re- 
quirement of any model is that the structural fluctuations 
should also be correctly reproduced. In this context, we 
have presented new QENS measurements of the diffusive 
dynamics of Ge on the length scale of ^3 A, which is 
most pertinent to the structural issues on the scale of 
the nearest neighbor distance. 

To understand this new QENS data, we have also per- 
formed ab initio simulations of liquid and amorphous Ge, 
the latter highlighting features of the liquid phase. These 
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simulations not only reproduce accurately the quasi- 
clastic broadening but also IXS data of acoustic phonons 
at small q, which concerns the structure on a length scale 
> 6 A. In addition, we have found very good agreement 
for the static structure factor, the pair correlation func- 
tion and the diffusion constant, the latter enabling an 
activation energy of 0.7 eV to be determined for diffu- 
sion processes, which is characteristic of a strong liquid. 

The success of AIMD is in contrast to classical MD sim- 
ulations that use two and three body potentials that, al- 
though they reproduce reasonably well the average struc- 
tural properties, do not give an accurate description of 
the dynamic a 17 i 38 . The QENS data and simulations also 
allow theoretical models, which relate the spectral line 
width to the static structure factor, to be tested. We 
have found that the second frequency moment given by 
Eq. ([6|) does not give a good description of the line width 
in the whole g-range whereas expression ([7]) describes this 
width well around the main peak in the static structure 
factor. 

Simulations also allow models to be investigated at the 
atomic level. Analysing the AIMD trajectories shows an 
average coordination number of 6-7 for the liquid and 
a broad range of bond angles with an average value of 
100 degrees. In contrast to this, for amorphous Ge, the 



coordination number is close to 4 and the average bond 
angle is closer to the value for a symmetric tetrahedron, 
even when the liquid density is imposed. In terms of 
simple models of the structure, the AIMD simulation is 
approximately consistent with the /3-tin hypothesis^. 

In conclusion, QENS data constitutes experimental ev- 
idence that can differentiate between different theoreti- 
cal and numerical models. The ab initio model presented 
here is in very good agreement with all new and existing, 
structural and dynamic experimental data. Further de- 
tails of the liquid structure are most likely to be obtained 
from more in-depth investigation of the MD trajectories. 
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